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Abstract. Raman light-scattering experiments in the antiferromagnetic phase of the 
Cu2Te2C'5(Bri_ :c CL c )2 compounds are analyzed in terms of a dimerized spin model for the 
tetrahedral Cu-clusters. It is shown that the longitudinal magnetic excitation in the pure Br 
system hybridizes with a localized singlet excitation due to the presence of a Dzyaloshinskii- 
Moriya anisotropy term. The drastic change of the magnetic scattering intensities observed 
when a proportion of Br is replaced by CI ions, is proposed to be caused by a change of the 
magnetic order parameter. Instead of being parallel/antiparallel with each other, the spins in 
the two pairs of spin-i order perpendicular to each other, when the composition x is larger 
than about 0.25. 



The compounds Cu2Te20sBr2 and Cu2Te20sCl2 are spin-i Cu-systems with coupled spin 
tetrahedra. The tetragonal P4 crystal structure of the two compounds and the susceptibility 
measurements are presented in ref. £Q . While the high-temperature susceptibility approaches 
that of uncoupled spins, at low temperatures the susceptibility is reduced and goes through a 
maximum at, respectively, 30 and 23 K. The maximum indicates that the 2 4 spin-^ states of 
the Cu-ions in the (distorted) tetrahedral clusters are dimerized so to create a singlet ground 
state separated from the excited states by a gap of about 40 K. Both the susceptibility and the 
heat capacity measurements clearly indicate a transition to an ordered phase at a temperature, 
11.4 K in the Br and 18.2 K in the CI system well below that of the maximum. Some of 
the bulk properties and the Raman-scattering results of the Br compound have been analyzed 
successfully by Gros et al. 0j in terms of a dimerized model for the four Cu spins of one 
tetrahedron, as determined by the following Hamiltonian: 

Ht = J\ (Si + S 2 ) • (S 3 + S 4 ) + J 2 (Si • S 2 + S 3 ■ S 4 ) (1) 

J x and J 2 are both positive. Defining the ratio r = J 2 / J\ then the energy differences between 
the singlet ground state Isj) and the excited states are: 2(1 — r)J 1 to a singlet |s 2 ), J l to a 
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one triplet state 1^), (2 — r)J 1 to two degenerate triplets \t 2 ) and |t 3 ), and 3J 1 to a quintuplet 
\q). This tetrahedral unit is coupled to the neighboring ones, and the Heisenberg interaction 
assumed in ref. 0] is 

< S F r) = - J o M z {SI + SI - St -St)+ 2 J c Ml , M 2 = \(Sf + SI - SI - SI) (2) 

within the mean-field (MF) approximation, i.e., the order parameter is one where the spins 
on each of the pairs 1-2 and 3-4 are parallel, but antiparallel with respect to the other pair 
of spins. This ordering takes full advantage of the ^-interaction on the expense of the J 2 - 
interaction. In the model of Gros et al. |2] for the Br compound r — 0.66 is smaller than 1, 
and the system is close to quantum criticality with a coupling parameter J c = 0.85^, only 
13% larger than the critical value Jc qc) — 0.75 J 1 . The model was investigated analytically 
neglecting the modifications due to higher lying levels. The high-temperature susceptibility 
is accounted for using J 1 = 47.7 K, the Neel temperature increases with an applied field, and 
the frequency of the longitudinal magnetic excitation in the ordered phase is estimated to be 
close to that observed by Raman spectroscopy. Here we shall present a complete MF analysis, 
which includes the calculation of the susceptibility and the bulk magnetization in the ordered 
phase and the field dependence of T N . It is performed numerically accounting for all effects 
of the total level scheme. The Raman-scattering cross section is derived within the RPA, and 
the field-induced |s 2 ) singlet excitation in the Raman spectrum is determined to be due to 
an intrinsic Dzyaloshinskii-Moriya (DM) anisotropy term in the Hamiltonian of the Cu spins 
and not to a DM contribution to the Raman-scattering operator as proposed in 0]. 

The bulk properties of Cu2Te20s(Bri_ :r Cl E )2 change gradually with x. This suggests 
that the only modification occurring is an increase of the effective J c in proportion to x, 
corresponding to the increase of T N from 11.4 to 18.2 K. However, the Raman scattering 
experiments, and to some extent also the magnetization measurements in the ordered phase, 
indicate that the nature of the ordered phase is substantially modified when x becomes larger 
than about 0.25. After a thorough analysis of alternatives to the M z order parameter in |3Jl, 
we conclude that the only one which is in accord with the experiments is the following 

^mf' = ~Jc V ^xy{S\ — S% — S% + S*!) + 2J^ y M xy , M xy = -^(Sf — S% — S% + SJ) (3) 

The M xy parameter describes an ordered state in which each pair of spins, 1-2 or 3-4, is in an 
antiparallel configuration, and the spins on the two different pairs are perpendicular to each 
other. When each pair of spins is antiferromagnetically ordered, the free energy contribution 
of the J 1 interaction cancels out to leading order. Thus the angle between the spins on the two 
different pairs is left undetermined by the intra-tetrahedral interactions in eq. whereas 
the DM-anisotropy favors the perpendicular configuration. 

The MF/RPA theory and the Raman- scattering cross section. - The ground state prop- 
erties and the excitations of the tetrahedral Cu-systems have been analyzed in terms of a 
general MF/RPA theory, see for instance [S]. The MF order parameter is determined in a 
self-consistent fashion, from a numerical diagonalization of the total MF Hamiltonian followed 
by a calculation of the thermal expectation value of the order parameter. The linear response 
at zero or non-zero frequency is derived from the non-interacting susceptibility tensor x°(w), 
with the components Xab( w ) determined in terms of the operators A and B, see section 3.3 in 
[5]. Introducing the MF-coupling constant J c in (2) then the 3x3 susceptibility tensor x(w), 
corresponding to A and B equal to the three components of S x + S 2 — S 3 — S 4 , is determined 
from the non-interacting one according to x( w ) = X°( w ) [ 1 — (J c /^)x • Introducing a 
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new, independent operator D, then the response Xdd( u ) * s derived from the 4x4 suscepti- 
bility tensor defined in the vector space consisting of the three Cartesian components plus D 
as the fourth component. Hence XzjdO^) IS ^ ne 44 component of the interacting susceptibility 
tensor determined from the non-interacting one as above, except that (J c /4) is replaced by 
a diagonal matrix with a zero in the 44 position and (J c /4) in the remaining part of the 
diagonal. 

The Raman scattering cross section R may according to Brcnig and Becker [5] be written 
T^im a hn( e i ' %m)( e o ' n i m )^i ' ^m! where and e c are the unit vectors of the incoming and 
outgoing electric field and n im are the unit vectors connecting exchange-coupled sites. When 
the incoming and outgoing light are polarized along the c axis, then R is proportional to 

]T [(Sl + S|) • (S| + S|) + a (S\ + S|) ■ (Sj-^ Sj" 1 ) + a (S 3 + Sj) • (S\ +1 + S 2 +1 ) + ■ • •" 

i 

R ^ A J2 [(Sl + S|) • (S| + Sl) - 2aM ■ (SI + S| - S 3 - S* 4 ) 

i 

This RPA expression includes all neighboring exchange-coupled pairs via the introduction of 
an effective a. Defining the operators D = (Si + S 2 ) • (S3 + S4) and A = S lz + S 2z — S 3z — S 4z 
then the scattering intensities in the Raman experiments are determined by 

T ^ K 1 _ e _L/fc B r Im {Xz»M + ^ 2m Ixaa{u) - 2aM z [ Xad {u) + XdaM] } ( 5 ) 

If a = J C /2J 1 then -R RP a commutes with the MF Hamiltonian. This corresponds to the 
case considered by Brenig and Becker of (3 — b. The Raman-scattering intensity vanishes 
identically in this case, independent of the value of M z . The value of a is unknown, and in 
the following we simply assume that Xdd( u ) 1S tne dominating term, i.e., that a ~ 0. 

Analysis of the experiments. The final model for Cu2Te205Bi - 2 includes J 11 J 2 , J c , 
and J* v defined by the equations ©-© plus two additional parameters. Jj is the coupling- 
parameter for the uniform magnetization component, 

W Z = -(J / M / + 2 A1B H).(S 1 +S 2 + S 3 + S 4 )+2J / M / 2 , M f = 1(8, + S 2 + S 3 + S 4 ) (6) 

of importance in the presence of an external field H. The S" 4 symmetry of the tetrahedra allows 
different kinds of anisotropic interactions. Based on the analysis of the Raman-scattering 
experiments we shall concentrate on the following DM anisotropy term: 

H A = D e [(S 1 -S 2 )x(S 3 -8 4 )] c (7) 

proportional to the c component of the cross product. The values of the five (six) parameters 
in the case of Cu 2 Te20 5 Br 2 are determined to be ( Jf y < 2.17^ =>• M xy = 0) 

J 1= 47.5K, J 2 = 0.7J 1; D C = 0.06J 1 , J c = 0.856J 1; ^ = 0.6^ (8) 

The strongest coupling in the system J 1 is derived from the high-temperature susceptibility 
measurements assuming fixed values for the remaining ones. The value of J j is mainly de- 
termined by the high-field magnetization measurements in the ordered phase, see fig. ^ but 
still has a large uncertainty of about 20%. Due to the D c anisotropy the bulk susceptibility 
becomes anisotropic. The anisotropy is minute in the paramagnetic, but substantial in the 
ordered phase. Independent of its sign D c implies that the ordered staggered moment M z in 
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Fig. 1 - High-field magnetization measurements after ref. [JJ. The smooth lines are the results of the 
model calculations. 



Fig. 2 - The bulk susceptibility at low temperatures. The open circles are the experimental results 
for polycrystalline samples [50. Single-crystal results with field parallel and perpendicular to the c 
axis are shown (solid symbols) for the case of the Br compound @J. All experimental results have 
been scaled so that they approach the theoretical high-temperature value of 0.7503/T[K] emu/mol. 
The solid/dashed lines are the calculated results when the field is parallel/perpendicular to the c axis. 




Fig. 3 - The Neel temperature as a function of applied field. The dashed lines (in the Br case) show 
the calculated T N when the field is perpendicular (upper one) or parallel (lower one) to the c axis. 
The solid lines are the averaged value to be compared with experimental results 2 4 obtained from 
polycrystalline samples (the solid squares) . 



eq. © is directed along the c axis. The bulk susceptibility along this direction is predicted to 
vanish exponentially at low temperatures whereas the perpendicular component stays approx- 
imately constant below T N . This is in qualitative agreement with experiments, as shown in 
fig. The experimental data are unfortunately somewhat disturbed by other contributions, 
as reflected in, for instance, a sample dependence of the results. The Zeeman energy gained 
in the perpendicular case induces a spin-flop transition, when the field is parallel to the c 
axis, to a phase where the staggered moments become perpendicular to the field. The model 
predicts the spin- flop field to be 3.7 T, nearly independent of the temperature below T N . This 
spin-flop transition has not yet been observed, but the anisotropy measurements in fig. [3 were 
made at a field of 5 T indicating that the spin-flop field is larger than 5 T. J c has no influence 
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Fig. 4 - Raman spectra of Cu2Te20sBr2 at 2.1 K in a magnetic field applied perpendicular to the c 
axis. The experimental results (symbols connected by dashed lines) I2HUII] have been obtained with 
two different resolutions and the theoretical intensities (solid lines) have been folded with Gaussians 
with the specified values of a. The calculated Raman shifts have all been subtracted by 1.6 cm -1 . 



on the bulk susceptibility, above T N , and is adjusted so that the model predicts T N = 11.4 
K at zero field. Figure shows the transition temperature for a polycrystalline sample as a 
function of field. It is remarkable that T N increases with the field. 

The Raman spectra of the C^T^C^Bri^CL;^ compounds have been studied in detail 
2,3,1]. The longitudinal magnetic excitation at zero wave vector is observed in the Br system 
at 2.1 K, when the incoming and scattered light are polarized along the c axis. This part of 
the results is shown in fig. El The high-resolution data revealed an additional peak at slightly 
higher energy than the longitudinal excitation, appearing when a field is applied perpendicular 
to the c axis. Gros et al. 0] interpreted the extra peak to be the \s 2 ) singlet level, which is 
becoming visible due to a DM component in the cross section combined with the level mixing 
induced by the field. The problem with this explanation is that the field-induced mixing is 
weak, and a DM cross section as large as the normal part would only produce a peak of the 
size of a hundredth of the main one. The alternative explanation is that it is the Hamiltonian 
itself, which is responsible for the mixing of the levels. The only possible term producing 
such an effect is the D c coupling defined by cq. J7J. D c gives rise to a coupling between the 
longitudinal mode and the |s 2 ) level even at zero field. The application of a field perpendicular 
to the c axis diminishes the energy difference between the two levels and the hybridization 
is strongly increased. The two peaks would attain equal intensities at a field slightly larger 
than the experimental maximum field. Figure 0] compares the intensities predicted by the 
model, eq. JHJ with a — 0, with the experimental ones. In these calculations D c and J 2 , 
plus one intensity scale parameter in each of the two cases, have been utilized as adjustable 
parameters. The experimental Raman shifts are subjected to an absolute uncertainty of the 
order of 4 cm -1 and a relative one of the order of 1 cm . The relative uncertainty has been 
reduced by performing a translation of each spectrum so that the phonon frequencies in each 
case coincide with those observed in the sample at zero field. Because of the large absolute 
uncertainty we have tried to mimic the relative behavior rather than to get the right values 
for the absolute frequencies. The model thus leads to a good account of both the frequencies 
and the intensities if shifting all the calculated peaks rigidly by —1.6 cm -1 . 
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Fig. 5 - Raman spectra obtained for small concentration x of CI ions at zero field and 4K, after ref. [7j. 
The arrows indicate the positions of magnetic scattering peaks. The intensity of the longitudinal 
magnetic excitation rapidly decreases with x and new peaks develop at higher frequencies. 

The model we shall present now for the Cii2Te205Cl2 system is more speculative. In 
comparison with the Br compound the a and c lattice parameters in the CI system are reduced 
by 3 and 1%, respectively, whereas the intra-tetrahedral distances between the Cu atoms are 
increased by about 1%, see ref. [Q. In correspondence to these changes J 1 is reduced in the 
CI compound whereas the Neel temperature is increased. A first guess of a model for the CI 
compound would be the same one as discussed above with some appropriate adjustments of the 
coupling parameters. The model based on this assumption predicts a strong Raman-scattering 
peak at about 40 cm -1 , two to three times as intense as in the Br system if assuming the 
scattering coefficients a lm to be the same. However, the Raman spectra clearly show a quick 
and systematic reduction of the magnetic scattering with the concentration x of the CI ions, 
see fig. One possible explanation would be that it is the Raman scattering circumstances, 
not the model itself, which is drastically changed. This could be achieved by introducing a 
non-zero a in eq. ijljl. which a should then approach J C /2J 1 in the pure CI system. We have 
analyzed this possibility, but found it to be both implausible and not capable of reproducing 
the observed non-monotonic changes of the Raman spectra with x. 

The analysis of the magnetization measurements shown in fig. ^indicates that r = J 2 / ' J 1 is 
close to 1 in the CI compound, so although J 1 is reduced J 2 increases with the CI concentration. 
When r ~ 1 then J c has to be about 2 J 1 in order to produce a transition to the M z ordered 
phase at T N = 18.2 K. Another possibility is that in both systems J c « J l whereas J£ v fss 
1.8J 1; then the increase of r from about 0.7 in the Br to about 1 in the CI compound is a 
sufficient cause for a change of the order parameter from M z to M xy . Based on this possibility, 
the model of Cu 2 Te20 5 Ci2 is proposed to be (J c < l.95J 1 => M z = 0): 



J 1= 40.7K, J 2 = J D D C = 0.06J 1: J% V = 1.8J V J* = 0.6^ (9) 



This model produces the magnetization and susceptibility results displayed in the figs. ^ 
and |3 and T N is found to be nearly independent of an external field, as shown in fig. 
Most importantly, the assumption of a M xy instead of a M z ordering implies a reduction of 
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the calculated magnetic Raman-scattering intensities by nearly a factor of 10. Furthermore, 
instead of one longitudinal mode the present model predicts two main peaks (at 43 and 75 
cm -1 ), plus a third one (at 57 cm -1 ) observable if the polarization vectors of the light have 
a component perpendicular to the c axis. This is, at least qualitatively, in much better 
agreement with the Raman experiments 0EI than that predicted by a M z ordered phase. 

Conclusion. - The analysis by Gros et al. 0] of the tetrahedral spin-^ cluster system 
Cu 2 Te 2 05Br 2 has been extended, and the Raman-scattering results have been explained quan- 
titatively by introducing a Dzyaloshinskii-Moriya anisotropy term in the Hamiltonian, which 
term also accounts (partly) for the anisotropy of the system below T N . The ordered phase 
in this compound is the one, where (S 1 + S 2 ) is parallel and (S 3 + S 4 ) antiparallel with the 
c axis. A spin-flop transition is predicted to occur when applying a field along the c axis. 
It is worth pointing out that the magnetic Raman spectrum (leaving out the two-magnon 
scattering part) depends on the resulting mean field, but not on the relative orientation of 
the magnetic moments on neighboring tetrahedra. A determination of this would require a 
neutron-diffraction investigation. The present analysis of the Raman spectra indicates that 
the ordering in Cu 2 Te 2 05Cl 2 compound is of the J 2 -type, where (S 1 — S 2 ) and (S 3 — S 4 ) are 
perpendicular to each other and to the c axis. The evidences are indirect and the J 2 -structure 
needs to be verified by diffraction experiments. 

The mean-field/random-phase approximation is a relevant starting point, whenever the 
interactions are strong enough to induce magnetic ordering. This is also true close to quantum 
criticality (Pr-metal may serve as a good example for that, jH)). The non-mean-field behavior of 
the present frustrated system appears to be accounted for by the treatment of each tetrahedral 
Cu-spin cluster as a local dimerized unity. The analysis shows that the specification of the 
properties of the local spin clusters requires at least three parameters J 1: J 2 , and D c . The 
values of the three MF parameters J*, J c , and J* y , which correspond to three different 
combinations of the exchange interactions between neighboring tetrahedra, indicate that both 
the interactions perpendicular and parallel to the c-axis chains of tetrahedra are important 
and therefore that the system is three- rather than one-dimensional. The good coincidence 
between the predictions of the MF/RPA theory and the experimental results is a further 
argument for a three-dimensionality of the magnetic interactions in the two compounds. 

* * * 
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